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Abstract — We study a two-tier macrocell/femtocell system 
where the macrocell base station is equipped with multiple 
antennas and makes use of multiuser MIMO (spatial multi- 
plexing), and the femtocells are "cognitive". In particular, we 
assume that the femtocells are aware of the locations of scheduled 
macrocell users on every time-frequency slot, so that they can 
make decisions on their transmission opportunities accordingly. 
Femtocell base stations are also equipped with multiple antennas. 
We propose a scheme where the macrocell downlink (macro- 
DL) is aligned with the femtocells uplink (femto-UL) and, Vice 
Versa, the macrocell uplink (macro-UL) is aligned with the 
femtocells downlink femto-DL). Using a simple "interference 
temperature" power control in the macro-DL/femto-UL direction, 
and exploiting uplink/downlink duality and the Yates, Foschini 
and Miljanic distributed power control algorithm in the macro- 
UL/femto-DL direction, we can achieve an extremely attractive 
macro/femto throughput tradeoff region in both directions. We 
investigate the impact of multiuser MIMO spatial multiplexing 
in the macrocell under the proposed scheme, and find that large 
gains are achievable by letting the macrocell schedule groups of 
co-located users, such that the number of femtocells affected by 
the interference temperature power constraint is small. 

I. Introduction 

It is widely recognized that spatial reuse is the single most 
valuable resource to dramatically increase the throughput of 
wireless cellular networks. However, deploying a very dense 
cellular infrastructure, with base station (BS) density that 
grows linearly with the user density, is not viable for a number 
of obvious practical and economical reasons. On the other 
hand, user-deployed WLANs (e.g., IEEE 802.11) achieve such 
dense spatial reuse in the unlicensed band. This solution has 
the advantage of providing very high data rates for short-range, 
mostly in-home, communication, but does not handle mobility 
as efficiently as cellular systems. Therefore, licensed cellular 
systems are naturally evolving towards two-tier architectures, 
where a large number of user-deployed femtocells operate 
under a common macrocell "umbrella", that fills in the gaps 
of the small cells tier and supports mobility. A large body of 
theoretical and standardization studies on this topic has been 
produced in recent years (for a small sample, see [1J, [2J, 0, 

a, 0). 

Most existing works focus on continuous transmission of 
macrocell user terminals (macro-UTs) and on the calculation 
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of the pdf of the signal-to-interference plus noise (SINR) at 
given receivers, where the tail of the SINR pdf is related to 
the probability of outage. This approach disregards the fact 
that the forthcoming generation of cellular systems (notable, 
3GPP LTE and IEEE 802.16m) is based on TDMA/OFDM, 
where dynamic scheduling is used in the macrocell tier for the 
the downlink (macro-DL) and for the uplink (macro-UL). For 
a macro-BS equipped with M transmit antennas and serving 
up to M macro-UTs on any time-frequency slot, the set of 
served macro-UTs (and therefore their location in the cell) 
may change on a slot by slot basis. This gives a statistical 
multiplexing opportunity to the femtocell tier: in the macro- 
DL slots, only the femtocells in the vicinity of a served macro- 
UT create significant interference to the macrocell tier; in 
the macro-UL slots, only the femtocells in the vicinity of an 
active macro-UT suffer from significant interference from the 
macrocell tier. 
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Fig. 1. Frame structure for the proposed cognitive femtocell system. 

In order to exploit the implicit statistical multiplexing due 
to the macrocell dynamic scheduling, in [6] we proposed a 
"cognitive" approach to femtocells, where we assume that 
the femto-BSs and the femto-UTs can decode the macro-BS 
allocation map for both UL and DL. Fig. [T] shows a possible 
arrangement of macro and femto frames allowing cognitive 
operations, in the case where both tiers use time division 
duplexing (TDD). Assuming that the positions of all terminals 
are known, [j] the femtocells can regulate their transmit power 
in order to guarantee a given "interference temperature" to the 
macro-UTs. In Section |Hl we review the details of a scheme 

'Femtocells are deployed in fixed positions, that can be made available 
through a database. Macro-UTs are mobile, but their position changes suffi- 
ciently slow so that through GPS and radio localization their position can be 
provided at a slow rate as protocol side information by the macro-BS itself. 



proposed in |6), based on linear beamforming and UL/DL du- 
ality, and we extend it to the case of multi-antenna macro-BS. 
Then, in Section|IIl]we discuss the coexistence of the multiuser 
MIMO (MU-MIMO) spatial multiplexing in the macrocell 
tier with the cognitive femtocells. In fact, it is intuitively 
clear that there is a tradeoff between the macrocell and the 
aggregate femtocell throughput: if the macrocell serves many 
macro-UT using spatial multiplexing, correspondingly many 
femtocells have to turn down their transmit power because 
of the interference temperature requirement, and therefore the 
femtocell throughput is decreased. In contrast, if the macro-BS 
serves only one macro-UT at each time-frequency slot, only a 
few femtocells are affected by the power control requirement 
but the macrocell tier does not exploit the full multiplexing 
gain and its throughput is decreased. We shall see that this 
problem is alleviated by scheduling approximately co-located 
groups of macro-UTs. 

II. SIMO/MISO INTERFERENCE CHANNEL 

We consider a single macro-BS with M antennas, serving 
K < M macro-UTs. In the same coverage area, a set of 
femtocells share the same frequency band. Both tiers operate 
in TDD. The channel gains are formed by two components: 
a pathloss factor constant in time (over a large number of 
time-frequency slots) and frequency-flat, and a time-frequency 
selective small-scale fading that changes independently on 
a slot by slot basis. For simplicity, we focus here on a 
single subcarrier. By symmetry of the fading distribution, our 
results extend directly to an OFDM system with independent 
scheduling on each subcarrier. |^] 

In the proposed scheme we have two types of slots: macro- 
DL/femto-UL and macro-UL/femto-DL (see Pig . |2]» . The fem- 
tocells operate in TDMA. Therefore, the number of femto- 
UTs actually present in each femtocell is irrelevant, since only 
one of them is active at any given slot and, for the sake 
of simplicity, it is sufficient to consider a single femto-UT 
per femtocell. Notice that the femtocells form a SIMO/MISO 
interference channel, coupled with the vector broadcast (DL) 
and multiaccess (UL) channel corresponding to the macrocell. 

Macro-DL/Femto-UL slot: Macro- and femto-UTs are 
equipped with a single antenna. The received signal at the 
fc-th macro-UT is given in general by 
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where h mc j 6 C ' J ^ is the channel vector from the macro-BS 
antenna array to macro-UT i, v, is the corresponding macro- 
BS beamforming vector, C denotes the set of all femtocells, 
hf.k is the scalar small-scale fading coefficient from femto- 
UT / to the macro-UT k, and z k ~ C7V(0, 1) is AWGN. The 

2 In general, a slight improvement can be obtained by scheduling jointly 
across the subcarriers. However, the improvement due to this more sophisti- 
cated multiuser scheduling across frequencies is marginal with respect to the 
system throughput achieved by the proposed scheme even with the suboptimal 
per-subcarrier scheduling. 
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Fig. 2. Interference channel model for the two-tier network with multiple 
antennas at the femto-BSs and macro-BS. The thick red lines indicate the links 
that determine the femto user transmit power, via the interference temperature 
power control. 



coefficients g(a, b) indicate pathloss between points a and b, 
as detailed in Section [III] The macro-BS is located at (the 
origin of the cell coordinate system). The macro-BS calculates 
the beamforming vectors as functions of the matrix 
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formed by the K active macro-UTs, enumerated without 
loss of generality from 1 to K. This can be obtained either 
by TDD reciprocity (open-loop) or by explicit channel state 
feedback [7|. In particular, here we consider Linear Zero- 
Forcing Beamforming (LZFB), such that v, is given by the 
i-th column of the Moore-Penrose pseudo-inverse of H mc 
normalized to have unit norm. Hence, h[^ c fe v; = for all 
i k. The macro-BS is subject to a total transmit power 
equal to Po, equally allocated over the K DL data symbols 

The femto-UT in the /-th femtocell transmits with power 
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Pf, regulated such that the interference caused at 



all active macro-UTs users is less than the target interference 
temperature k. Hence, we have 
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where Pi is the peak femtocell power. The SINR for macro- 
UT k is given by 



SINR™ 



-DL 



g(M)|h£c. fc v*l 2 £ 

-E /eC <?(^/)l^/l 2 ^/ 



(4) 



and the corresponding instantaneous rate for macro-UT k on 
the current slot is given by R k = log(l + SINR" 1C ~ DL ). For 
simplicity, we assume that the macro-BS schedules at each slot 
K out of U K macro-UTs, picked at random with equal 
probability. Hence, by averaging over the fading realization 
and the ( K ) sets of macro-UTs, we obtain the sum-throughput 
of the macrocell tier in the DL. When K = M and the fading 
is Rayleigh i.i.d., using the results in [8], this can be given in 
closed form for fixed U macro-UT positions. 



The femto-BS receivers have L antennas each, and make 
use of linear MMSE detection, which maximizes SINR over 
all linear receivers. The received signal vector at a femto-BS 
/ is given by 
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The linear MMSE receive vector for estimating the desired 
symbol Xf from yf is given by u/ = otfUJ hfj where a/ > 
is chosen such that ||u/|| = 1, and X/ is the interference- 
plus-noise covariance matrix in |5]), given by 
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The receiver forms the scalar observation yf = u^y/, and the 
corresponding SINR is given by 



SINR^ UL = P f h i }j'EJ 1 h 



/,/• 



(7) 



Similarly to what argued before, the instantaneous rate of 
femtocell / is given by R f = log(l + SINR^- UL ). By 
summing over all the femtocells and averaging over the fading 
and the (^) active macro-UTs sets (notice that they have an 
influence through the femtocell transmit powers Pf), we obtain 
the sum-throughput of the femtocell tier in the UL. 

Macro-UL/Femto-DL: On the macro-UL/femto-DL slot, 
each femtocell has multiple antennas whereas the macro-UTs 
have single antenna each. Insisting on linear beamforming 
strategies, each femto-BS sends the L-dimensional signal vec- 
tor xy = WfSf, where denotes the transmit beamforming 
vector and Sf is the corresponding (coded) data symbol for its 
own intended femto-UT. 

The received signal at the macro-BS is given by 
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The BS forms the scalar observation y k = r^y^ for detecting 
Xfe, where is the receive beamforming vector for macro-UT 
k. 

The received signal at the femto user in femtocell / is given 

by 
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Calculating the instantaneous SINRs SINR™ C_UL and 
SINR^ DL fr om &8t and from d9b, respectively, is straightfor- 
ward. In particular, from UL/DL duality (see |9|), extended 
to the MISO/SIMO interference channel in 1151 . we have that 
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where Q mc .k denotes the transmit power of macro-UT k and 
Q f denotes the transmit power of the femto-BS /. The power 
allocation across the macro-UTs and the femto-BSs depends 
on the realization of the path losses and small scale fading 
components (which determine the beamforming vectors). 

A. Implementation issues 

In order to calculate the beamforming vectors u/, using 
the matrix inversion lemma we can write u/ = fif'K.~j 1 h.fj, 
where j3f > is another normalizing proportionality constant 
and K/ is the received signal covariance matrix given by 

K f = H f +g(f,f)h f , f hl f P f (11) 

Hence, the MMSE beamforming vectors can be conveniently 
calculated by using a sample covariance estimate of K/, 
from the whole received femto-UL slot, and an estimate of 
the desired signal channel hfj obtained by using UL pilots 
symbols, as in standard coherent detection for MIMO channels 
(e.g., currently implemented in IEEE 802.1 In). 

The other practical implementation problem of the proposed 
scheme consists of calculating the transmit powers Q mc ,fc 
and Qf in the macro-UL/femto-DL slot, for fixed unit-norm 
beamforming vectors w/ = Uf and = v^. We propose to 
use the well-known Yates-Foschini-Miljanic distributed power 
allocation algorithm (see ifTTI . lfT2l ). that is guaranteed to 
converge to the solution. 

For all femtocells /, fix the target DL SINR 
7^ C - DL = SINR^ UL . For all active users k fix 
the target UL SINR 7j: nc " UL = SINR^ C " DL . Let 
SINR^ DL ({Q / },{Q mCjfe },{u / }) denote the femtocell 
/ DL SINR for fixed beamforming vectors and transmit 
powers, and let SINRr~ UL ({Q/}, {Qmc.fe}, {u/}, {v fc }) 
denote the macro-UT k UL SINR for fixed beamforming 
vectors and transmit powers. Then, the iterative distributed 
power control algorithm, in our case, is given by: 

1) Initialization: let n = and let Q^l k = Pq/K, — 
P f for k = 1, . . . , K and all / e C' 

2) Iterations: for n = 1,2,3,... do 
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In order to implement this scheme, a sequence of adjacent slots 
should be allocated to the same group of K macro-UTs, and 
at each slot the receivers measure their SINR and report their 
measurements to the transmitters such that the power values 



by letting Wf = Uf and r-fc — it is possible to achieve can be updated according to (12i 



TABLE I 
Simulation Parameters 



Parameter 


Notation 


Value 


Macro Cell Side Length 


L 


1000 m 


Path Loss Parameter 


8 


50 m 


FC Radius 


ro 


10 m 


Distance between two FC 


I 


40 m 


Path Loss Exponent 


a 


3.5 


Wall Partition Loss 


w 


5 dB 


Min SNR at cell edge 


SNR min 


10 dB 


Number of antennas at macro-BS 


M 


8 


Number of antennas at femto-BS 


L 


5 



III. Numerical results 

In line with [6|, we consider a unit-side square cell 
[-1/2,1/2] x [-1/2,1/2], with the macro-BS located at 
the origin 0, and F 2 femtocells are centered at points of 
coordinates ( 2t ~^ +1 , 2j '^g +1 ), for i,j = 0,...,F - 1. 
Femtocells are disk-shaped with radius rf c , shielded from the 
outdoor environment by walls. For two points a, b, the distance 
dependent path loss component is given by 



g(a, b) 
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where d(a, b) denotes the distance between a and b modulo 
the centered unit square (torus topology); n(a, b) counts the 
number of walls between points a and b (i.e., n(a, b) = if 
both a and b are outdoor or they are in the same femtocell, 
n(a,b) = 1 if either a or b is indoor (inside a femtocell), and 
n(a,b) = 2 if a and b are in different femtocells); w is the 
wall absorption factor; 5 is the"3 dB" pathloss distance; a is 
the outdoor pathloss exponent. 

We fix the macro-BS power Pq such that the received SNR 
(without interference) for a macro-UT at the cell edge is 10 
dB. By varying the value of the interference power temperature 
k and letting Pi = 30 dB^jwe obtain the Pareto boundary of 
the throughput tradeoff region achievable with the proposed 
scheme in the macro-DL/femto-UL slot, as described before. 
The tradeoff region for the dual channel, i.e., macro UL/femto 
DL is obtained by using the same beamforming vectors and the 
iterative power control algorithm. We distinguish between the 
cases of colocated and non-colocated macro-UTs. In the first 
case, we macro-BS schedules K users roughly located in the 
same position of the cell, such that they are separated enough 
to have independent small-scale fading, but they have the same 
pathloss with respect to the macro-BS and all the femtocells. In 
the second case, the K macro-UT positions are independently 
selected with uniform probability over the cell. Fig. [5] shows 
the comparison of the Pareto boundaries of the tradeoff regions 
for the colocated and non-colocated case (supremizing over 
K), showing a clear advantage for the colocated case, wich 
is possible when the macro-UT density is large enough so 
that K approximately colocated users can be found. In all 

3 Assuming a peak rate in a femtocell in isolation is 10 bit/s/Hz, corre- 
sponding to uncoded 1024 QAM, we have SNR = 101og 10 (2 10 - 1) = 
10*log 10 (1023) « 30dB. 
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Fig. 3. Throughput tradeoff region comparison for colocated users in macro- 
DL/femto-UL. The different colors indicate the tradeoff region for different 
K. 
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Fig. 4. Throughput tradeoff region comparison for non-colocated users in 
macro-DL/femto-UL. The different colors indicate the tradeoff region for 
different K. 



the results, femtocells are assumed to be "open access", 
therefore, macro-UTs are located only outdoor since a macro- 
UT inside a femtocell would be automatically "swallowed" 
by the femtocell, and served as a femtocell user. We already 
showed in (6| that the impact of closed-access femtocells 
is minimal, in contrast to what observed for conventional 
"legacy" systems, thanks to the proposed cognitive scheme 
and the interference temperature power control. 

Figs. [3] and |4] show the throughput tradeoff region (femtocell 
sum throughput vs. macrocell sum throughput) achieved by the 
proposed scheme in the macro-DL/femto-UL slot (averaged 
over random user positions), for colocated and non-colocated 
macro-UTs, respectively. As the number of served macro- 
UTs increases, the macrocell throughput increases initially 
up to a certain maximum value (in our case, the highest 
macrocell throughput is obtained for K = 6 users) and then 
decreases. This can be expected from the typical behavior 
of linear LZFB precoding. For the non-colocated case, the 
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Fig. 5. Throughput tradeoff region comparison (Pareto boundary) between 
colocated and non colocated users in macro-DL/femto-UL. 
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Fig. 6. Throughput tradeoff region for K = 6 colocated users in macro- 
UL/femto-DL for increasing iterations of the power control algorithm. 

femtocell throughput decreases significantly as K increases, 
since more and more femtocells are affected by the power 
control interference temperature limitation. This is because 
K random macro-UTs positions are selected at each slot. 
Instead, for colocated macro-UTs, the value of K has no 
effect on the femtocell throughput, therefore, K — 6 achieves 
(approximately) uniformly best performance over the whole 
throughput range. 

Fig. [6] and shows the throughput tradeoff region for the 
macro-UL/femto-DL slot, for K = 6 and different iterations of 
the power control algorithm. In every iteration of the power 
control algorithm, the peak power constraint is imposed for 
both the femtocells as well as the macro user terminals^] 
We notice that already for 3 iterations the algorithm yields 
almost the target "dual" macro-DL/femto-UL region, and for 
6 iterations the region is indistinguishable. Notice also that 
because of the imposed peak power constraint, duality does 

4 In the macro-DL, each user is given an equal fraction of the power Pq . In 
the macro-UL, each user is constrained to use a power no greater than Pg / K. 



not strictly hold. Nevertheless, as seen from these plots, the 
impact of the peak power constraint on the macro-UL/femto- 
DL slot is basically negligible for this realistic range of system 
parameters. 

IV. Conclusions 

Overall, the ergodic rate region achievable by the proposed 
scheme is very competitive with other schemes proposed or 
analyzed in the current literature, considering that it can be 
achieved with a very simple protocol and low-complexity 
signal processing. For example, operating the system at the 
achievable throughput tradeoff point with macrocell through- 
put of 15 bit/s/Hz and femtocell throughput of 1000 bit/s/Hz, 
we can achieve 600 Mb/s over 40 MHz of system bandwidth 
of average macrocell symmetric data rate (both UL and DL), 
and 64 Mb/s per femtocell over the same system bandwidth (in 
our system geometry we have 625 femtocells per macrocell). 
Given these rather outstanding numbers, we believe that the 
proposed scheme is an attractive option for "beyond 4G" 
future wireless networks. Of course, several issues need further 
investigation, as for example the system operations and per- 
formance with multiple macrocells, the protocol overhead for 
implementing cognitive femtocells and adapting the power by 
the iterative algorithm, and the effect of scheduling co-located 
macrocell users on the macrocell user channel correlation, 
which may limit the effective macrocell multiplexing gain. 

References 

[1] S. Rangan, "Femto-Macro Cellular Interference Control with Sub- 
band Scheduling and Interference Cancelation," Arxiv preprint 
arXiv: 1007.0507, 2010. 

[2] V. Chandrasekhar, J. Andrews, and A. Gatherer, "Femtocell networks: 
a survey," IEEE Commun. Mag., vol. 46, no. 9, pp. 59-67, 2008. 

[3] S. Yeh, S. Talwar, S. Lee, and H. Kim, "WiMAX femtocells: a perspec- 
tive on network architecture, capacity, and coverage," IEEE Commun. 
Mag., vol. 46, no. 10, pp. 58-65, 2008. 

[4] P. Xia, V. Chandrasekhar, and J. Andrews, "Open vs closed access 
femtocells in the uplink," Arxiv preprint arXiv: 1002.2964, 2010. 

[5] H. Dhillon, R. Ganti, F. Baccelli, and J. Andrews, "Modeling and 
analysis of k-tier downlink heterogeneous cellular networks," Arxiv 
preprint arXiv:1103.2177, 2011. 

[6] A. Adhikary, V. Ntranos, and G. Caire, "Cognitive femtocells: Breaking 
the spatial reuse barrier of cellular systems," in Information Theory and 
Applications Workshop (1TA), 2011. IEEE, pp. 1-10. 

[7] G. Caire, N. Jindal, M. Kobayashi, and N. Ravindran, "Multiuser mimo 
achievable rates with downlink training and channel state feedback," 
Information Theory, IEEE Transactions on, vol. 56, no. 6, pp. 2845- 
2866, 2010. 

[8] H. Kim, S. Lee, K. Lee, and I. Lee, "Transmission schemes based 

on sum rate analysis in distributed antenna systems," Arxiv preprint 

arXiv: 1105.6199, 2011. 
[9] P. Viswanath and D. Tse, "Sum capacity of the vector Gaussian broadcast 

channel and uplink-downlink duality," IEEE Trans. Inform. Theory, 

vol. 49, no. 8, pp. 1912-1921, Aug. 2003. 
[10] F. Negro, I. Ghauri, and D. Slock, "Beamforming for the underlay cog- 
nitive MISO interference channel via UL-DL duality," in 5th Int. Conf. 

on Cognitive Radio Oriented Wireless Networks and Communications 

(CROWNCOM). IEEE, pp. 1-5. 
[11] R. Yates, "A framework for uplink power control in cellular radio 

systems," Selected Areas in Communications, IEEE Journal on, vol. 13, 

no. 7, pp. 1341-1347, 1995. 
[12] G. Foschini and Z. Miljanic, "A simple distributed autonomous power 

control algorithm and its convergence," IEEE Trans. Veh. Tech. , vol. 42, 

no. 4, pp. 641-646, 1993. 



